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Abstract 

The energy spectra of spin-1/2 electrons under two-dimensional magnetic 
field modulations are calculated beyond the one-band approximation. Our 
formulation is generally applicable to a modulation field with a rectangu- 
lar lattice symmetry. The field distribution within a plaquette is otherwise 
arbitrary. The spectra being obtained are qualitatively different from their 
electric modulated counterparts. Peculiar features of the spectra are that, 
for an electron with a g factor precisely being equal to two, no matter how 
strong the modulation is, the zero-energy level seems to be unaffected by the 
modulation and is separated from higher energy levels with a nonzero energy 
gap. Moreover, there is a two-fold degenerancy for all states with positive 
energies with respect to spin flip. These features agree with earlier analytical 
studies of the periodically magnetic modulated systems. 

PACS numbers: 71.25-s, 73.20Dx, 74.60.-w 
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I. INTRODUCTION 



Due to the progress of the submicron technology, one begins to observe quantum be- 
haviors in the transport measurements of the field modulated two-dimensional electron gas 
(2DEG). Much effort has been devoted to the study of the spectral and the transport prop- 
erties of a 2DEG in a periodic magnetic field with a nonzero uniform background Bq. These 
studies can be divided into two classes, depending on whether the field modulation is one 



For a one-dimensional magnetic modulated (1DMM) system, there are two characteristic 



length scales: the magnetic length, A = ya/e£>o, associated with the uniform background 
Bq and the period a of modulation. By varying the ratio A/a, the electron mobility and 
the magnetoresistance oscillate between extrema.S The oscillating behavior of the latter is 
similar to the Weiss oscillation in a one-dimensional electric modulated (IDEM) 2DEG.0 
These oscillations manifest the variation of the band widths: the (longitudinal) conductivity 
is proportional to the width of the Landau-level (LL) broadening, which is an oscillating 
function of both A/a and LL indices due to the field modulation. Under special conditions 
(the so-called flat band conditions), the band width can be zero and electrons become 
immobile. Besides the transport property, other aspects of the 1DMM system have also 
been studied, such as collective excitations,! inelastic light scatterings,! surface states^ 
effect of electron-electron interactions,! and the effect due to an additional two-dimensional 
electric modulation.! The 1DMM systems may also be used as spin polarizers for magnetic 
dipoles.i Recently, the 1DMM systems have been realized experimentally by covering a 
regular array of superconductor^ or micromagnet! on the top of a 2DEG, in which the 
observed magnetoresistance oscillation agrees very well with the theoretical prediction. 

For the two-dimensional magnetic modulated (2DMM) systems, the Landau levels are 
not only broadened but also split to several subbands with an intricate fractal structure.00 
This is similar to the Hofstadter spectra for the two-dimensional electric modulated (2DEM) 
systems.0 In fact, within the one-band approximation, both spectra are exactly the same 
to linear order in modulation fields.00 For a square lattice in the one-band approximation, 
a Landau band is split to p subbands when there are p/q (p and q are relative prime 
integers) flux quanta per plaquette. Previous calculations of the 2DMM systems are either 
restricted to the one-band approximationjll or to the multiband calculation but with a 
specific magnetic field such that the flux quantum per plaquette is one (or one half). Bill 
Other aspects of the 2DMM systems have been studied, such as the collective excitations^ 
and the degeneracy of the ground states.01l! There are some recent attempts to observe the 
peculiar transport properties due to the fractal band structure. However, as far as we know, 
this goal has not been achieved for the 2DMM systems.0 

In the present work, the energy spectra of a spin-1/2 2DEG under two-dimensional 
magnetic field modulations are calculated beyond the one-band approximation, in which 
the Zeeman term is also included. In most of the above studies, the Zeeman effect is 
not included, however. For the electric modulated systems, neglect of the Zeeman term 
is justified, because the periodic electric field does not couple to electron spin and this 
term only contributes to a constant energy shift, (g e eh/4m)Boa z , where g e is the electron 
g factor and a z is +1(— 1) for spin- up (spin-down) electrons. However, this is not the case 
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for the magnetic modulated systems. After deriving the multi-band Harper equation, which 
is generally applicable to magnetic modulations with arbitrary strength and shape, as long 
as the field has a rectangular lattice symmetry, we show that the inclusion of the Zeeman 
term leads to qualitative changes in the energy spectra. Particularly, when g e = 2, the 
most disparate result occurs for the lowest energy level — it is not broadened by the field 
modulation, and is separated from higher energy bands by a finite gap. Moreover, there 
exists a two-fold degenerancy for all states with positive energies with respect to spin flip. 
These results agree with earlier mathematical analysis of the 2DMM systems. 

The paper is organized as follows: the multi-band formalism is presented in Sec. II; the 
band structure is presented in Sec. Ill; and Sec. IV is devoted to summary and discussions. 



II. MULTI-BAND FORMALISM 

A. Magnetic translation symmetry 

We consider a 2DEG under the influence of a magnetic modulation with a rectangular 
symmetry. The Hamiltonian is 

a(r)) + 9 jB(r)o- z , (1) 

where A (r) and a(r) are the vector potentials for the uniform background field B and the 
modulation field 6(r) = B(r) — B , respectively. In this paper, unless specified explicitly, 




we choose A = ^h/eBo as the unit of length, %u c as the unit of energy (lo c = eBo/m is the 
cyclotron frequency), and B as the unit of magnetic field. In the absence of modulation, 
the Hamiltonian Hq can be solved exactly with eigenvalues = n + 1/2 + g e a z /Ax^ H can 
be expanded as H + Hi + H 2 + H a , where Hi and H 2 are the terms linear and quadratic 
in the vector potential a(r), respectively, and H a = (g e /4)b(r)a z is the modulated Zeeman 
term. The vector potential a(r) can be Fourier decomposed as a(r) = 2 g ^o a g e * gr ; where 
g are the reciprocal lattice vectors of the rectangular lattice. By choosing the Coulomb 
gauge, the Fourier components a g are equivalent to ib s g x z/g 2 , where b s are the Fourier 
components of b(r) and g — |g| (not to be confused with the electron g factor g e ). It is 

convenient to rewrite the exponential e ig r as e* g '^e* g ' R , where the electron coordinate r is 
decomposed into a fast-moving cyclotron coordinate £ and a slow-moving guiding-center 
coordinate R = r — £. (See Appendix |A].) Then it can be shown that 

ff i = -E3Vl^ E - (2) 



g 2 d\ 



Similarly, 
and 



jEEg-g'^^ 4 ^ 11 . (3) 
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Ha = ^J2 & g e^e ig - R (4) 

Due to the underlying magnetic translation symmetry^ of the Hamiltonian, it is con- 
venient to diagonalize the Hamiltonian on a basis which respects this symmetry. The 
unperturbative basis can be constructed as follows. By choosing a Landau gauge with 
Ao(r) = (— y, 0), the magnetic translation operators are 

rp^ _ e aid/dx ^ rp^ _ e a 2 (d/dy-ix) ^ /g\ 

where a\ and a 2 are the lattice constants for the rectangular lattice. It is not difficult to 
show that, if there are p/q flux quanta per plaquette with an area a\a 2) then T\, T$, and 
H mutually commute. This is also true for the unmodulated Hamiltonian H Q , of course. 
Thus we can construct the explicit form of the magnetic Bloch states for H Q , which are the 
common eigenstates of H , T% and T 2 9 :E1 

°° 2-7T 

|n,k)= E die-^^l^h 1), (6) 

l=-oo fl l 

where n is the LL index, k = (k\, k 2 ) is the magnetic Bloch momentum, d\ are complex 
coefficients that are periodic in I with period p (i.e., di +p = di), and \n, k\) are the common 
eigenstates of H and T\. 

Since both H and Hq are diagonal with respect to k, it is clear that for the modulation 
part H' = H — H , we have (n, k.\H'\n', k') = (n, k\H'\n', k)5k,k'- Therefore, the a-th 
eigenstate of H can be written as 

OO OO c\ 

l«,k) = E E <?e- <(9/p)w |n,*i - -/>, (7) 

n=0/=-oo fll 

where the unknown coefficients d^j, as di above, are periodic in I with period p. Basically, 
the strategy below is to diagonalize the Hamiltonian matrix on the unperturbed basis and 
to solve for its eigenvalues. 

In deriving the matrix elements of the Hamiltonians in Eqs. (2)-(4) on the unperturbed 
basis, the expression (n, k% — 2nl / ai\e ls '^e ls ' K \n' , ki — 2ixl'/a\) will be encountered frequently, 
thus we focus on its derivation first. Firstly we rewrite the exponential in a slightly different 
form, e ls '^e 1 ^, to connect with the magnetic translation symmetry, where g = g x z 
and R = R x z. Since the two dynamical variables, £ and R, decouple and operate on 
different parts of the Hilbert space (see Appendix the matrix elements of e lg '^e 1 ^ can 
be evaluated with the help of Eq. ( |Al| ) . The result is 



n, ki I 



2-7T 

n'M V) = 5 hl ,_iP kll (g)U nn/ (g), (8) 



where g = (gi,g 2 ) = {2Ttl/a l ,2nm/a 2 ), P kl i{s) = e - 7TiI,:h ^ p e 27Tikl ^ a2 e- Ml ^ p ', and U n 



(n 



e ts '^\n'). The magnetic flux condition, a\a 2 = 2irp/q, has been used. 
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B. Multi-band Harper equation 



When the energy eigenstate |a, k) is expanded on the basis of \n,ki — 2nl/a-\) using 
Eq. (|7p, the eigenvalue equation, H\a, k) = E a (k)\a, k), takes the following form, 



ESPd® + E e^'W-Hn, h - — \H'\n', h - = f^M 

ai 



n',/' 



(9) 



where I = pr + s and /' = pr' + s'. (r, s, r', and s' are all integers such that < s, s' < p.) 
Firstly we need to calculate the matrix elements of H'. With the help of Eqs. (@)-(fD an d 
(0), we have 



n, 


h- 


2tt/ 


Hi 


i 

n , 
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= Si,i'-i 
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n. 
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2nl'\ 


= 4 6l > 








ai / 





dU nn ,(g\ ) 
' OX 



A=l, 



(10) 



Combining Eqs. (§) and C[T0|) , we finally obtain 



n' g^O 
1 



l dU nn >(g\) g e cr ZTT , , 

U=l H tVnri'lgJ 



, fen- &cr' 



d (Q) 



(a) 



- o E E E s • g'i^(g + g')^(g + g')< J 
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where P ks (g) = e- 2mfc2r/ai P felS (gl = e ik -te x2 ) e - i3l92/2 e- i92(2 ™ /ai) . It is a multi-band gener 

tiorE 



alization of the Harper equationllH (see Eq. (0)), thus Eq. is called as the multi-band 
Harper equation.il It has to be solved in conjunction with the following identities concering 
inter-LL transitions (for n > n'), 



dU nn >(g X) 



A=l 



n\ \^/2 

9- 

m \y/2 



(„ _ „' _ g 2 /2)L«7 n ' - g 2 L n n T^ +1 



(12) 



where g>_ = g± — ig 2 and L~"^ n {UL\ = 0) are the associated Laguerre polymials with the 
argument g 2 /2. 

Until now, no approximation has been used. This multi-band Harper equation applies 
to general shape of magnetic field distribution with a rectangular lattice symmetry. To 
simplify the calculation, from now on, we assume the spatial field modulation is 6(r) = 
26 10 [cos(27ra;/a) + cos(27n//a)] with a square lattice symmetry (i.e., a x = a 2 = a). In the 
numerical calculations, the eigenvalues are obtained by diagonalizing a (1 + n cnt )p x (1 + 
n cut )p matrix, where the cut-off n cut has to be large enough to ensure the eigenvalues being 
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obtained converge to the correct result. The precise spectra that include the effect of inter-LL 
transitions are shown in the next section. 

Before closing this subsection, we show that, under the so-called one-band approximation, 
Eq. flllD can indeed be reduced to the usual Harper equation. When the inter-LL transitions 
and the terms quadratic in b s are neglected, a^j — > did" (due to the periodicity, di +p = d~i, 

and then Eq. ([0]) is reduced to 



there are only p independent coefficients, i.e., do 
the one-band equation: 

M n (i^j |4_ ie 2 " fc2/a + d s+1 e' 2mk2/a + 2d s cos 



, dp. 





£l 0) 



(13) 



where M n (q/p) = (6 10 /2)[L^ + L} 1 _ 1 + (g e a z /2)L n }e~ 9 ^o/ 4 is an overall factor that scales the 
energy, and g\ = (2-n/a) 2 = 2nq/p. Apart from the factor M n (q/p) where the spin-related 
term is included, Eq. (|T3|) is precisely the same as the Harper equation for a 2DEM system.0 
The spectrum for E n (q/p)/M n (q/p) within the one-band approximation is thus trivial:0lll 
irrespective of the LL index n, it is the usual Hofstadter spectrum calculated for a 2DEM 
system. There is one exception, however. When n — 0, g e — 2, and a z = — 1, M (q/p) is 
equal to zero, and then E (q/p) - 
all. Actually, the equality E (q/p) 



E^, as if the field modulation exerts no influence at 
E^ is valid even beyond the one-band approximation. 



This is discussed in more details in Sec. Ill 



III. FRACTAL BAND STRUCTURE 

In this section we show the band structures for both weak and strong modulations. The 
influence of the Zeeman term is particularly emphasized. For an unmodulated 2DEG with 
no Zeeman effect, the energy spectrum consists of discrete, dispersionless LLs.Eil These LLs 
are highly degenerate because of the continuous translation symmetry, which gives an infi- 
nite degenerancy, and the spin-flip symmetry, which gives a two-fold degenerancy. When a 
periodic modulation is introduced, which breaks the continuous translation symmetry, we 
expect that the degenerancy for each LL will be lifted. Indeed, in the one-band approxi- 
mation, one finds that each LL is broadened and split to several intricate energy subbands. 
The way these subbands split is the same for every LL in the one-band approximation (see 
the discussion at the end of the last section). However, when the inter-LL transitions are 
included, the exact calculations shown below reveal that the subband structures are actually 
different for different LLs, thus lead to much more complicated structures. 

In Fig. 1, the spectrum of electrons with g e = under a weak square modulation field with 
frio° 2 — O.20 o (in the usual units, where (fi = h/e is the flux quantum) is shown. Because 
there is no Zeeman splitting, there is no need to distinguish the spin-up electrons from the 
spin-down electrons and the spectrum for only one spin direction is shown. The calculation 
is done with a cut-off n cut = 9. The result with a larger cut-off at n cut = 14 shows no visible 
difference from Fig. 1. Notice that the abscissa is the inverse of the magnetic flux, q/p. In this 
and the following figures, it is assumed that, while changing the magnetic flux by varying B , 
the modulation amplitude bio is fixed. For the weak modulation case, the envelope for each 
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energy band is largely determined by the scaling factor M n (q/p). Obviously, some features 
specific to the one-band approximation no longer exist. For example, the interband couplings 
remove the symmetry of the butterfly diagram. Similar effect of symmetry breaking is also 
observed in the multi-band calculation for the 2DEM systems.il 

The two-fold degenerancy for electron spin is lifted when the Zeeman effect is not negli- 
gible. In our 2DMM systems, the Zeeman effect does not only give an energy shift, but also 
induce inter-LL transitions (see Eq. flTU| ) for H a ). Thus, the interplay between the orbital 
effect (Hi + H 2 ) and the Zeeman effect (H a ) leads to different spectral structures between the 
spin- up and the spin-down electrons, which are shown in Fig. 2. The spectrum in Fig. 2(a) 
(Fig. 2(b)) is for a spin-down (spin- up) electron with a g factor equals to one. It can be seen 
that, in addition to the overall constant Zeeman energy shift g e a z /A due to the background 
field, the spectra show qualitative differences from Fig. 1. For example, the structures for 
the second lowest energy bands near q/p = 0.9 in Figs. 2(a) and (b) are visibly different 
from that in Fig. 1, and are different from each other. 

The spectrum in Fig. 3(a) (Fig. 3(b)) is for a spin-down (spin-up) electron with a g factor 
equals to two. For this particular g factor, two significant features are observed. The first is 
that the lowest energy level for a spin-down electron is flat and equals to zero to very high 
precision. The second is that the positive-energy spectrum is degenerate with respect to 
spin flip, whereas the flat band in Fig. 3(a) has no counterpart in Fig. 3(b). Both features 
persist for stronger modulations. 

The two-fold degenerancy for the positive-energy states in the case of g e = 2 indicates 
that there must be an additional symmetry even in the presence of the Zeeman term. It 
was pointed out by Aharonov and Casherc^l that this degenerancy results from a symmetry 
transformation which simultaneously changes the direction of the electron spin and the spa- 
tial dependence of the wavefunction. This symmetry can be related to the supersymmetryii 
or to the chiral (75I invariance by connecting our problem to the (1 + l)-dimensional theory 
of Dirac fermionsl^l 

Moreover, based on an abstract mathematical analysis, Dubrovin and Novikov showed 
that, for a spin-down electron with g e = 2, there always exist the zero-energy states in 
the 2DMM systems, no matter how strong the modulation is.0 Furthermore, by using 
topological arguements, they proved that, although the continuous translation symmetry 
is broken in the 2DMM system, the degenerancy of this zero-energy states is the same as 
that for the unmodulated system. One may wonder whether this unexpected degeneracy is 
symmetry- related and, if it is, what is the nature of this symmetry. In fact, it was shown by 
Gendenshtein that, because the Hamiltonian in Eq. (1) with a e = 2 can be factorized into 
a product of two conjugate first-order differential operators,^ such a symmetry is indeed 
present for the zero-energy states.@ However, this is not a symmetry of the original Hamil- 
tonian, but rather of one of the two first-order differential operators. (For more details, see 
Ref. [0].) 

In comparison with the g e ^2 cases, there is one more unique feature for the g e = 2 case, 
which becomes obvious when the modulation is quite strong. Figs. 4(a) and (b) show the 
strongly modulated band structures of spin-down electrons with g e = 1 and 2, respectively. 
The modulation strength is b\oa 2 = O.80o- For the case of g e = 2, besides the fact that the 
lowest energy level remains flat despite the strong mixing between the unperturbated LLs, 
it is apparent that the zero-energy level is isolated from the intermingled fractal structure 



7 



with a finite energy gap. This is true for even stronger modulations. On the contrary, such a 
behavior does not appear in Fig. 4(a) for the case of g e = 1. A simple explanation of the gap 
above the flat band is as follows:0 if the gap collapses at a particular modulation, such that 
a state from a higher energy band merges with the zero-energy one, then the degenerancy 
of the zero-energy states will increase by one. However, this is impossible because this 
degenerancy in the modulated system must be the same as that in the unmodulated one, as 
mentioned above.0 Therefore, the flat band has to be separated from higher bands. 

Fig. 5 shows the dependence of the energy bands on the g factor when the total flux per 
plaquette is B a 2 = ltfio and b w a 2 = O.20 o - In this figure, the sign of the g factor refers 
to the direction of spins (+ for up, — for down). (Note that the sign convention for g e 
applies only to Fig. 5, but not to the previous discussions.) It is clear that the lowest band 
is broadened as soon as g e ^ —2. Notice that the spectrum for g e = 2 and that for g e = —2 
are identical, except for the lacking of the flat band at zero energy for g e = 2. In addition 
to the zero-energy flat band, the band widths of other energy bands can also shrink to zero 
at some particular values of g e . For example, The width of the second lowest energy band 
is zero when g e = —0.12. However, unlike the shrinking of the lowest band at g e = —2, this 
'pinch' point moves if a different flux value of B^a 2 is chosen. 

IV. SUMMARY AND DISCUSSIONS 

In this paper, we present an accurate multi-band calculation of the energy spectra of the 
2DMM spin-1/2 electronic systems, in which the Zeeman effect is also taken into account. 
We find that, when the Zeeman energy is not negligible, the spectra are changed qualitatively 
with respect to their electric modulated counterparts. Moreover, in the special case when 
the electron g factor is two, it is found that: 1. the positive-energy eigenstates have a two- 
fold degeneracy with respect to spin flip; 2. for the spin-down electrons, the ground states 
seem to be unaffected by the periodic modulation and remains highly degenerate even in 
very strong modulations; 3. the ground states are separated from higher energy states with 
a finite energy gap. However, these special properties no longer exist if the g factor is not 
equal to two. 

In real systems, the electron g factor can be changed by either varying the width of 
a quantum well that holds the electrons,^ or by applying a hydrostatic pressure^ to the 
sample. Thus, it is possible to combine these two methods and to design an experiment in 
which the g factor can be continuously tuned around two.@ Under such a circumstance, the 
peculiar and robust spectral properties of the lowest energy level should exhibit itself through 
the transport properties. Since the band splitting may suppress the band conductivity^ it 
is advised to keep the flux per plaquette at one or a simple fraction in order to observe the 
conductivity enhancement induced by level broadening away from g e = 2, or, conversely, the 
conductivity reduction at g e = 2. 

However, most of current experiments can achieve only weak modulations. For example, 
consider a source of a periodic magnetic field, B(r) = B -\-(B /2)(cos2iTx/a+cos2ny/a) > 0. 
The field modulation felt by a 2DEG at a distance d below the source becomes Bd(r) = 
B + (B e~ 27Td/a /2)(cos2nx / a + cos27n//a).0 This corresponds to b w = e - 27rc( / a /4 in our 
calculation. For typical values such as a = 1 /im and d = 10 nm, 6 10 is equal to 0.23. The 
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optimum value of 610 is 1/4, when d = 0. It can be larger than 1/4 only if the amplitude 
of the modulation field is larger than the background field. In this case, the total field B(r) 
reverses direction in some regions. 

In the future, to bring theoretical results much closer to the real experiments, ingredients 
such as disorder and electron-electron interaction have to be included in the calculation. It 
may also be necessary to include an extra two-dimensional electric modulation, which is 
inevitably induced due to the strain exerted by the ferromagnetic or the superconductor 
grid at low temperatures in recent experiments.&uiiil It will be very interesting to investigate 
the influence of these factors on the energy spectra reported here. 

ACKNOWLEDGMENTS 

This work is supported by the National Science Council of Taiwan under contract No. 
87-2112-M-007-008. 

APPENDIX A: GUIDING-CENTER COORDINATE 

In the semiclassical calculations of the transport properties of the quantum Hall systems, 
in which the applied magnetic field is strong and the disorder potential is smooth,!! one 
usually decomposes the electron coordinate r into a fast-moving cyclotron coordinate, $, = 
z x (p + A ), and a slow-moving guiding center, R = r — £. In the present study, although 
the periodic field variation is not required to vary smoothly, the derivation of the multi-band 
Harper equation can be simplified with the help of this decomposition. 

It can be shown that (^1,^2) and (Ri,R 2 ) are independent conjugate pairs respectively, 
i.e., [£1,^2] = ~h [-^1,-^2] = i, and Rj] = for i,j = 1,2. Thus the exponential e 8gr can 
be decomposed as e* g '^e* s ' R = e ig '£e* s R , where g = gxz and R = Rxz. For a Landau gauge 
with Ao(r) = (— y,0), we have £ = (—id y ,id x + y) and R = (— id x , —id y — x). Therefore, 
the magnetic translation operators in Eq. (5) can be rewritten as Tj = e iR i a i, j = 1,2. 
Consequently, we have the following very useful identities, 

Ti \n, kx) = e iRiai \n,h) = e ikiai \n,h), 

T 2 \n,kx) =e ik2a2 \n,ki) = \n,k x -a 2 ). (Al) 

The second equation is a direct result of the commutation relation between Ri and R 2 , 
[Ri,R 2 ] = i. Therefore, the 'rotated' guiding-center coordinate is the generator of the 
magnetic translation. It can be verified that e tRiai \n, k) = e tklCll \n, k) and e tR2qa2 \n,k) = 
e iqk2a2 \n,k) using Eq. (|A1|) and the periodicity of d~i (see Eq. @). Thus, \n, k) are indeed 
the common eigenstates of H , T 1; and T 2 . Similarly, due to the periodicity of d^j (see 
Eq. ([?[)), one can prove that the energy eigenstates \a, k) are also the eigenstates of T\ and 

rpq 

J 2 • 
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FIGURES 



FIG. 1. The lowest five energy bands for spin-down electrons with g e = under a weak square 
modulation field. The modulation strength is fixed at &ioa 2 = O.20o while the total flux per 
plaquette, Boa 2 , varies. The energy is in units of ?ko c and the inverse flux per plaquette is in units 

of^o 1 - 

FIG. 2. The lowest five energy bands for a spin-down (a) and a spin-up (b) electrons with 
g e = 1. The modulation strength and the units being used are the same as those in Fig. 1. 

FIG. 3. The lowest five energy bands for a spin-down (a) and a spin-up (b) electrons with 
g e = 2. The modulation strength and the units being used are the same as those in Fig. 1. It 
can be seen that the lowest energy level for a spin-down electron is flat. Also, the spectrum is 
degenerate with respect to spin flip, except for the zero-energy flat band in (a). 

FIG. 4. The lowest five energy bands for a spin-down electron with g e = 1 (a) and 2 (b). The 
modulation strength is &ioa 2 = 0.8<^o- 

FIG. 5. The dependence of the band widths on the g factor when the total flux per plaquette 
Boa 2 = 10o and b\oa 2 = 0.2^o. The sign of the g factor refers to the direction of spins (+ for up, 
— for down). The energy bands at g e = —2 are indicated by vertical bold lines. 
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